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Growth  Through  Agrlcuttural  Progress 


ESTIMATING  DEPENDABLE  ANNUAL  STREAMF^W 
IN  THE  UNGLACIATED  ALLEGHENY  PLATEAU  i/ 


D.  L.  Brakensiek 


2/ 


The  BI5  Problem  Area  in  Soil  Conservation  can  be  viewed  essentially 
as  the  Appalachian  Plateau  Physiographic  Province.  Fenneman  {^)  z/   indicates 
that  this  Province  for  systematic  study  shoiild  be  divided  into  seven  areas. 
These  differ  in  the  character  of  their  underlying  rocks,  in  altitude,  in 
degree  of  dissection,  and  in  the  presence  or  absence  of  glaciation. 
Fenneman 's  division  was  partly  verified  by  Morisawa  {k)   through  the  use  of 
quantitative  morphological  parameters.   It  was  fo\ind  that  the  Allegheny 
Mountain,  Unglaciated  Allegheny  Plateau,  and  Cumberland  Plateau  Divisions 
showed  significant  differences  in  certain  morphological  parameters. 

The  study  reported  here  is  limited  to  the  largest  and  most  typical 
division  of  the  province,  the  Unglaciated  Allegheny  Plateau  Division.  Water- 
sheds utilized  in  this  study  further  define  the  area  of  study  and  future 
application  of  results  (see  Appendix  Figures  1  to  k).      Geographically  the 
study  area  covers  southeastern  Ohio,  southwestern  Pennsylvania,  western 
West  Virginia,  and  a  portion  of  eastern  Kentucky. 

Morisawa  (k)   describes  the  Unglaciated  Allegheny  Plateau  Division  as 
mature  in  stage,  with  fine  texture,  and  moderate  to  great  relief.   Both 
surface  elevation  and  relief  are  lowest  on  the  western  edge  in  Ohio  and 
Kentucky  where  elevation  averages  1,200  to  1,300  feet  with  relief  of  200  to 
300  feet.  Eastward  surface  elevations  reach  1,500  to  1,600  feet.   In  West 
Virginia  valleys  are  steep-walled  and  deep,  with  an  average  relief  of 
1,000  feet.   The  siirface  of  the  division  is  characteristically  rugged,  being 
mainly  in  valley  side  slopes.  Divides  are  long  and  linear  with  many  inter- 
mittent areas  of  rolling  upland.  Valley  flats  are  generally  narrow  and  con- 
fined to  larger  streams  and  lower  major  tributaries.   Stream  valleys  are 
generally  narrow,  steep-sided  and  often  deeply  incised.  However,  many  become 
gentle  in  gradient  near  their  heads  with  broadly  rounded  divides. 

BASIC  DATA 

Annual  streamflow  amo\ints,  in  vuiits  of  inches  over  the  watershed,  for 
watersheds  ranging  in  size  from  under  100  acres  to  i^-00  square  miles  were 
tabulated  on  a  calendar  year  basis.  For  consistency  the  calendar  year  was 

1/  Contribution  from  the  U.  S.  Hydrograph  Laboratory,  Soil  and  Water 

Conservation  Research  Division,  Agriciiltural  Research  Service, 
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2/  Hydraulic  Engineer,  U.  S.  Hydrograph  Laboratory,  Soil  and  Water 

Conservation  Research  Division,  Agricultural  Research  Service, 

U.  S.  Department  of  Agriculture,  Beltsville,  Md. 
3/  Numbers  in  parentheses  refer  to  Literature  Cited  at  end  of  this  report. 


used  with  tabulated  U.  S.  Weather  Bureau  annual  precipitation  data.  Agricul- 
tural Research  Service  watersheds  at  Coshocton,  Ohio,  and  U.  S.  Geological 
Survey  gaged  watersheds  were  utilized.  A  list  of  watersheds  is  given  in 
Column  1  of  Table  1. 

According  to  Climate  of  the  States  (9)  long-time  mean  annual  precipi- 
tation in  the  study  area  probably  varies  from  36  to  ^8  inches.  However,  since 
tabulations  of  annual  streamflow  data  were  often  for  different  periods  of 
time,  the  geographic  precipitation  difference  wo\ild  be  further  accentuated. 
Thus  to  correct  for  this  annual  watershed  precipitation  was  computed  for  use 
as  a  concomitant  factor. 

Annual  precipitation  amounts  in  inches  for  appropriate  stations  were 
tabulated  from  the  U.  S.  Weather  Bureau's  "Climatic  Summary  of  the  United 
States  --  Supplement  for  I93I  through  I952"  (7).  For  later  years  the 
annual  summaries  of  State  Climatological  Data  (8)  were  utilized. 

Iso-lines  of  annual  precipitation  were  constructed  covering  the  pertinent 
sections  of  each  state  for  each  year  of  record.  From  the  iso-line  maps 
annual  values  of  watershed  precipitation  were  determined  for  each  watershed, 
and  thus  from  these  average  annual  watersheds,  precipitation  was  calculated 
for  each  watershed's  length  of  record.   This  procedure  was  thought 
appropriate  because  very  seldom  was  there  one  or  more  raingages  actually  within 
a  watershed.   Column  3  of   Table  1  indicates  that  calculated  watershed  precipi- 
tation varied  from  37  to  57  inches.   This  variation  reflects  geographic  as 
well  as  length  of  record  differences. 

Obtaining  detailed  and  consistent  watershed  physiographic  parameters 
was  not  possible  because  different  watersheds  were  covered  by  different 
scaled  topographic  maps,  l/62,500;  l/2i<-,000j  or  l/4,800.  However,  the 
following  simple  items  were  obtained:   watershed  size  in  square  miles, 
perimeter  in  miles,  length  of  principle  water  course  in  miles,  and  maximum 
relief  in  feet.   More  detailed  physiographic  factors  would  have  required  a 
field  survey. 

The  exact  location  of  the  gage  site  for  each  watershed.  Column  1, 
Table  1,  can  be  found  in  the  U.  S.  Geological  Survey  Water  Supply  Papers  for 
the  Ohio  River  Basin. 

FORMULATION  OF  STUDY 

The  purpose  of  this  study  is  to  develop  a  procedure  for  estimating 
annual  streamflow  for  any  level  of  dependability.  Assuming  that  annual 
streamflow  amounts  are  statistical  variates,  we  can  write  a  formal  expression 
(1)  for  estimating  probability  -  magnitudes  as 

(^  =  Q  (1  +  CV^  K)  (1) 

where  for  a  particular  watershed 

Q^^  =  annual  streamflow  that  will  not  be  equalled  for  a  return 

period  T 
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Q   =  average  annual  streamflov 

CV  =  parameter  reflecting  inherent  Yariability  of  annual  streamflow 

K  =  standardized  deviate  for  a  probability  level  of  (i/t);  its 
value  also  depends  on  the  statistical  distribution  of  the 
variates  annual  streamflov 

The  application  of  equation  (l)  to  a  particular  watershed  requires  that 
three  items  be  known  or  procedures  for  estimating  them  be  available.   The 
items  are: 

(1)  An  appropriate  statistical  distribution  for  annual  streamflow. 

(2)  Q,  average  annual  streamflow. 

(5)  CVq  ,    inherent  variability  parameter. 

These  will  be  considered  in  the  next  sections. 

Statistical  Distribution 

No  clear-cut  statement  can  be  made  about  the  statistical  distribution 
most  appropriate  for  annual  streamflow  amounts.   The  choice  in  many  hydro- 
logic  studies  is  often  made  by  trying  several  and  selecting  the  one  that 
most  consistently  fits  the  available  sets  of  data.   This  study  is  no  excep- 
tion in  that  the  normal  distribution,  log-normal  distribution,  and  gamma 
distribution  were  all  investigated. 

Annual  streamflow  data  were  plotted  on  normal  and  log-normal 
probability  paper  to  test  whether  one  or  the  other  woiild  be  appropriate.  A 
graphical  test  of  the  adequacy  of  the  gamma  distribution  is  not  a  simple  task 
as  the  distribution  cannot  be  standardized  for  construction  of  a  "universal" 
probability  paper.  However,  several  test  graphs  were  prepared  for  this 
distribution. 

The  plot  of  each  set  of  watershed  data  on  normal  probability  graph  paper 
indicated  that  this  statistical  distribution  consistently  described  the  data. 
(See  Appendix  Figures  5  to  ^l).  A  discussion  of  the  graphical  scheme  is 
presented  in  Appendix  1. 

There  is  probably  some  theoretical  justification  for  a  normal  distribu- 
tion describing  the  accumulated  probability  fimction  of  annual  streamflow 
in  this  study  area.   One  csm  consider  that  annual  streamflow  is  composed  of 
the  sum  of  many  independent  streamflow  events  occurring  during  the  year. 
According  to  the  central  limit  theorem  then  this  sum,  annual  streamflow, 
under  very  general  restrictions  (2)  would  be  asymptotically  normally  dis- 
tributed. For  an  area  of  low  annual  streamflow,  where  very  few  independent 
events  occur  in  a  year,  annual  streamflow  would  be  skewed  considerably  from 
a  normal  distribution. 


Table  1.- 

-Watershed  Information. 

Length 

NAME 

q" 

p" 

^^^ 

CVp 

Size 

of 
record 

OHIO 

Ins. 

Ins. 

Sq^.  mi. 

Yrs. 

Short  Creek 

ik.ks 

4l.O 

0.3356 

0.1606 

124.0 

15 

Lisbon  Creek 

15.57 

58.6 

— 

— 

6.08 

10 

Mill  Creek 

13.61 

57.79 

.3997 

.1476 

27.5 

20 

Salt  Creek 

15.^0 

59.4 

— 

— 

10.6 

10 

Tar  Hollow  Creek 

11.60 

40.6 

— 

— 

1.50 

10 

Home  Creek 

11.15 

57.2 

.4208 

.1574 

1.64 

20 

Yellow  Creek 

15.42 

59.8 

.3262 

.1445 

148.0 

16 

Coshocton,  No.  97 

12.3 

57.4 

.5785 

.1474 

7.16 

20 

Coshocton,  No.  95 

12.0 

57.5 

.5854 

.1512 

4.02 

19 

Coshocton,  No.  9^ 

12.2 

57.5 

.3734 

.1512 

2.57 

19 

Coshocton,  No.  92 

12.3 

57.5 

.3782 

.1512 

1.44 

19 

Coshocton,  No.  5 

11.37 

57.5 

.5975 

.1512 

0.545 

19 

Coshocton,  No.  10 

9.61 

57.5 

.4o4o 

.1512 

.191 

19 

Coshocton,  No.  I96 

1U.65 

57.8 

.5801 

.1552 

.475 

20 

Coshocton,  No.  I85 

10.57 

56.5 

.4517 

.1648 

.116 

19 

Coshocton,  No.  177 

7.67 

56.5 

.4670 

.1561 

.118 

18 

Coshocton,  No.  172 

11.58 

56.5 

.4558 

.1561 

.068 

18 

WEST  VIRGINIA 

Buffalo  Creek 

20.6 

46.8 

.2988 

.l4l9 

115.0 

24 

West  Fork  River 

22.38 

48.1 

•  5246 

.1582 

181.0 

26 

Skin  Creek 

22.46 

49.0 

— 

— 

25.7 

11 

Buckhannon  River 

28.6 

57.2 

.2587 

.1425 

277.0 

26 

West  Fork  River 

25.5 

49.1 

— 

— 

102.0 

10 

Wheeling  Creek 

15.76 

41.2 

.5911 

.1551 

282.0 

16 

Little  Kanawha 

River 

24.38 

51.4 

.5146 

.1450 

155-0 

19 

Do 

21.7 

48.1 

.5257 

.1595 

586.0 

26 

Coal  River 

17.16 

45.5 

.2884 

.1421 

595.0 

26 

Little  Coal  River 

16.93 

44.5 

.5198 

.1452 

270.0 

26 

Steer  Creek 

17.63 

45.8 

.5506 

.1548 

166.0 

19 

^fe£Ldow  River 

24.91 

48.2 

.2154 

.1567 

287.0 

26 

PENNSYLVANIA 

• 

Conno  quene  s  s  ing 

Creek 

18.12 

59.0 

.2867 

.1401 

556.0 

25 

Crooked  Creek 

20.80 

44.5 

.2990 

.1458 

191.0 

19 

Dunkard  Creek 

17.14 

45.4 

.5509 

.1501 

229.0 

16 

South  Fork  Ten  Mile 

15.62 

41.16 

.5556 

.1566 

180.0 

25 

KENTUCKY 

Tygarts  Creek 

17.05 

42.59 

.4546 

.1675 

241.0 

17 

Johns  Creek 

16. 5 

44.4 

.5749 

.1654 

55.7 

15 

Triplett  Creek 

21.9 

45.2 

.5800 

.1662 

47.5 

16 

Poor  Fork 

25.1 

48.4 

.5006 

.1240 

82.1 

16 

Little  Sandy 

16.4 

45.6 

.4400 

.1989 

598 

19 

Levis a  Fork 

15.1 

45.5 

.5564 

.1598 

386 

18 

Cumberland  River 

25.7 

49.2 

.2820 

.1202 

374 

16 

Blaine  Creek 

15.2 

45.2 

.4600 

.21 22 

217 

18 

Licking  River 

15.7 

45.5 

.5997 

.1814 

l4o 

19 

Theoretically  the  gamma  distribution  with  zero  skew  woiild  also  have 
fitted  the  data.  H.  C.  S.  Thorn  (5^  6)  gives  meLximum  likelihood  estimators 
for  the  gamma  distribution.  The  data  of  Appendix  Figure  6,   when  substituted 
in  these  estimating  equations,  gave  a  skewness  statistic  of  about  O.76.   This 
is  not  zero,  but  if  the  confidence  limits  were  known,  they  might  encompass 
zero.  However,  the  point  is  that  for  this  set  of  data  the  gamma  distribution 
did  not  reduce  to  the  normal  distribution  which  seemingly  fits  the  data 
quite  well.  Even  if  it  had,  the  simplicity  of  the  normal  distribution  would 
make  it  preferable  in  this  study  to  the  gamma  distribution. 

Average  streamflow,  Q 

In  an  earlier  section  it  was  mentioned  that  average  annual  watershed 
precipitation  varied  from  watershed  to  watershed  because  of  geographical 
location  and/or  different  lengths  of  record.   Since  part  of  the  variation  of 
average  annual  streamflow  from  watershed  to  watershed  is  certainly  due  to 
the  precipitation  difference,  average  annual  watershed  precipitation  is  used 
as  a  concomitant  factor. 

SeveraJ.  watershed  physiographic  factors  were  also  studied  as  possible 
independent  variables  in  an  estimation  equation  for  Q.  Three  factors  were 
investigated:   watershed  size,  perimeter,  and  relief  ratio. 

It  was  found  that  watershed  size  or  perimeter  in  combination  with 
average  watershed  precipitation  were  significant  factors.   The  relief  ratio 
was  not  found  to  be  a  significant  factor.   Since  size  is  a  more  common  measiire- 
ment  than  perimeter  it  was  selected  for  the  final  estimating  equation.   The 
predicting  equation  is  as  follows: 

_  1.915    0.026 
Q  =  0.0133^+  (P)      (A)  (2) 

where 

Q  =  average  annual  streamflow,  inches 

P  =  average  annual  watershed  precipitation,  inches 

A  =  watershed  size  in  square  miles 

The  coefficient  of  determination  associated  with  the  logarithmic  transforma- 
tion of  (2)  is  R^  =  0.86.  The  standard  deviation  for  each  exponent, 

1.915  +  0.035 
and  0.026  +  0.001 

indicates  that  the  factors  are  quite  highly  significant  in  the  statistical 
sense .  In  a  later  example  computation  it  will  be  shown  that  the  seemingly 
small  exponent  of  the  area  term  is  also  of  considerable  practical  significance. 


Streamflov  variability,  CVg 

The  variability  factor,  CVq,  can  be  estimated  for  each  watershed  by 
applying  least  squares  procedures  to  equation  (l),  i.  e. 

(4)  =     ^(^/-^)^ 


X  k2 

This  is  possible  because  the  K-values  are  easily  determined  from  a  Table  of 
Areas  for  the  normal  distribution.  Accumulated  probability  values  are 
determined  from  the  plotting  position  formula 

m  /  (n  +  l) 

where 

m  =  order  number  after  items  are  ranked,  large  to  small 

n  =  length  of  record 

Equation  (l)  was  only  fitted  to  sets  of  data  of  at  least  15  years  of  length. 

It  was  felt  that  reliable  accvimvilated  probability  functions  could  not  be 

established  for  shorter  records.   The  fitted  lines  are  also  indicated  in 
Appendix  Figiires  5  "to  kl. 

The  calculated  CVq  values  were  related  to  easily  obtained  factors. 
The  best  and  most  simple  predicting  equation  found  is  as  follows: 

CVq  =  0.2540  -  0.0071  (Q)  +  1.5750  (CVp)  (5) 

where 

CVq  =  parameter  in  equation  (l) 

Q  =  average  annual  streamflow,  inches 

CVp  =  coefficient  of  variation  associated  with  an  appropriate  record 
of  annual  precipitation. 

The  coefficient,  of  determination  associated  with  equation  (3)  is 
r2  =  0.71.   Both  Q,  ard  CVp  are  at  least  highly  significant  factors.   CVp 
would  be  calculated  from  an  appropriate  annual  watershed  precipitation 
record  by  the  usual  formula, 

y7:(p  -  p)2/  (N  -  1) 

CVp  =   


^(P)/N 


In  terms  of  uncorrected  sums  of  squares  it  could  be  calculated  by  the 
following  equation: 


Q  co\ild  be  estimated  by  equation  (2). 


DISCUSSION 

Equations  (2)  and  (3)  together  with  the  assumption  of  normality  com- 
pletely define  the  probability-magnitude  estimation  equation  for  annual 
streamflow,  equation  (l). 

Equation  (2)  indicates  that  average  annual  streamflow  is  almost  related 
to  the  square  of  annual  precipitation.  Even  though  the  exponent  of  the  area 
factor  is  small  it  reflects  a  significant  trend,  that  is,  average  streamflow 
increases  with  increasing  watershed  size. 

Equation  (3)  points  out  some  interesting  items  about  year  to  year  vari- 
ability of  annual  streamflow.   Streams  that  have  higher  average  annual  flows 
exhibit  less  year  to  year  variations  in  flow.  Also,  watersheds  that 
experience  less  variation  in  annual  precipitation  also  have  less  variation 
in  annual  streamflow. 

The  following  niomerical  example  illustrates  the  application  of  the 
results  derived  in  this  study: 

GIVEN: 

(a)  Watershed  area   A  =  356  square  miles 

("b)  Average  annual  watershed  precipitation   P  =  39  inches 
(see  Table  2). 


(c)  CV  -     o.l4  (calculation  shown  in  Table  2). 


P 
COMPUTATIONS: 

(a)  From  equation  (2) 

9  =  0.0153^  (59)1-915  (356)0-026 

"q  =  17.25  inches 

If  the  area  term,  raised  to  the  0.026  power,  had  been  neglected  Q  would  have 
been  l4.8l  inches.  This  wo\ild  have  been  an  underestimate  of  over  I5  billion 
gallons  of  water.' 
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(b)  From  equation  (3) 

CVp  =  0.23^+0  -  (0.0071)  (17.25)  +  1.5750  (O.li^) 

CVp  =  0.3320 

Thus,  the  estimating  equation  for  dependable  annual  streamflow  for  this 
example  watershed  is 

^  =  17.25  (1  +  0.33K) 

The  above  equation  can  be  used  in  one  of  several  ways. 


Table  2. --Precipitation  Record  for  Connoquenessing  Creek, 
Hazen,  Pa. 


YEAR 

Pi/ 

COMPUTATIONS 

In. 

1932 

35 

33 

3^ 

>:  P  = 

9lh 

3h 

32 

35 

39 

r  p2  = 

3Q63h 

36 

ho 

37 

kS 

N   = 

25 

38 

37 

39 

36 

19^+0 

k2 

p  = 

^  p/n  = 

974/25 

i+l 

36 

k2 
h3 

kk 

p"  = 

39.0 

35 

37 

2 

kk 

j        2k 

25(38634)-(974)^ 

he 

46 
32 

CVp  = 

t'       (974)^ 

hi 

36 

kQ 

hi 

k9 

33 

CVp  = 

(1.02)   X 

V    (0.018103) 

1950 

k6 

51 

kl 

CVp  = 

0.14 

52 

k2 

53 

33 

5^ 

kh 

55 

3h 

1956 

51 

1/  Calculated  from  annual  iso-precipitation  maps 


Table  3. — Selected  Return  Periods  and  Associated 
Frequency  Factors. 


Prob.  of  a 

Retvirn 

Frequency- 

lesser  value 

period,  T 

factor,  K 

Percent 

Year 

95 

1.05 

1.6if 

90 

1.11 

1.28 

85 

1.18 

l.Oi^ 

80 

1.25 

.Qk 

75 

1.33 

.67 

70 

1.^3 

•  52 

65 

1.5i^ 

.59 

60 

1.67 

.25 

55 

1.82 

.15 

50 

2 

0 

i^5 

2.22 

-  .15 

ko 

2.50 

-  .25 

55 

2.86 

-  .59 

30 

3-33 

-  .52 

25 

k 

-  .67 

20 

5 

-  .Qk 

15 

6.67 

-l.oi<- 

10 

10 

-1.28 

5 

20 

-1.64 

2 

50 

-2.05 

1 

100 

-2.35 

Using  the  values  of  K  found  in  Table  5^  values  of  axmual  watershed  streamflow 
can  be  calculated  for  selected  return  periods.   On  the  other  hand,  it  may  be 
preferable  to  plot  the  equation  itself  on  normal  probability  paper  and  read 
values  for  any  return  period.  Figure  k2   shows  the  above  derived  equation 
plotted  on  normal  probability  paper.   Two  points  needed  to  orient  this 
line  can  be  calculated  as  follows: 

(1)  One  point  should  always  be  taken  at  Probability  =  50^  since  for 
this  value  K  =  0  and  (Q/Q)  =  1. 

(2)  The  second  point  could  be  calculated  for  an  accumulated  probability 
value  of  95  percent.   From  Table  2  the  value  of  K  for  this  probability  level 
is  K  =  1.64.   Thus,  the  ratio 

(095^)  /"q)  =  1  +  0.33  (1.64)  =  1.52 
The  second  point  then  has  coordinates 

(95^  1.52). 
Of  course,  any  value  read  from  a  graph  such  as  Appendix  Figure  42  wovild  be 
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multiplied  "by  the  mean  to  obtain  the  actual  Q  value.   In  the  example,  the 
values  read  from  the  ordinate  scale  would  be  multiplied  by  Q  =  I7.25. 

For  a  specific  application  it  might  be  feasible  to  read  an  average 
annual  watershed  precipitation  value  from  Appendix  Figures  1  to  4.  However, 
actual  precipitation  station  records  shovild  be  utilized  for  calculating 
the  coefficient  of  variation  of  watershed  annual  precipitation. 
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APPENDIX 

Any  graphical  test  for  the  "goodness  of  fit"  of  a  theoretical  accumulated 
probability  function  to  observed  data  reduces  to  a  plotting  of  the  observed 
reduced  variate  to  the  theoretical  reduced  variate.   That  is  to  say  the 
quantity 

(x  -  .U   )/cT 

is  plotted  versus  the  quantity 

(y  -  y)/  Sy 

Where  the  former  is  the  theoretical  reduced  variate,  in  this  study  read 

from  a  table  of  areas  of  the  normal  distribution,  and  the  latter  the 

observed  reduced  variate.   The  condition  of  equivalence  leads  to  the  following 

(y  -  y)/  Sy  =  (x  -  u  )/<r 

which  reduces  to  y  /  y  =  1  +  (s  /  y)(x  -  i^  )/c^ 

or 

(y  /  y)  =  1  +  (cv)K. 

This  is  the  same  as  equation  (l). 
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Mean  Annual  Precipitation,  Inches 


OHIO 


Based  on  period  1931-55 

laolinea  axe  drawn  tliroiigli  points  of  approximately  equ2d  value.    Caution  should  be  used 
in  interpolating  on  these  maps,  particularly  in  mountainous  areas. 

Figvire  1. 
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